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ABSTRACT 

Relaxed, massive galactic objects have been identified at redshifts z = 4,5, 
and 6 in hydrodynamical simulations run in a large cosmological volume. This 
allowed us to analyze the assembly patterns of the high mass end of the galaxy 
distribution at these high zs, by focusing on their structural and dynamical prop- 
erties. Our simulations indicate that massive objects at high redshift already 
follow certain scaling relations. These relations define virial planes at the halo 
scale, whereas at the galactic scale they define intrinsic dynamical planes that 
are, however, tilted relative to the virial plane. Therefore, we predict that massive 
galaxies must lie on fundamental planes from their formation. 

We briefly discuss the physical origin of the tilt in terms the physical processes 
underlying massive galaxy formation at high z, in the context of a two-phase 
galaxy formation scenario. Specifically, we have found that it lies on the different 
behavior of the gravitationally heated gas as compared with cold gas previously 
involved in caustic formation, and the mass dependence of the energy available 
to heat the gas. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: Evolution 

— galaxies: Formation — galaxies: Fundamental Parameters — hydrodynamics 

— methods: numerical 



^Grupo dc Astrofisica, Departamento de Fi'sica Teorica, Modulo C-15, Universidad Autonoma de Madrid, 
Cantoblanco E-28049, Spain 

^Department of Physics and Astronomy, University of California at Irvine, Irvine, CA 92697, USA 

■^Depto. de Fisica y A.C., Universidad Miguel Hernandez, E-03202 Elche, Alicante, Spain 



-2- 



Introduction 



One of the outstanding yet most important problems in astrophysics is how and when 
galaxies formed within the framework of the expanding Universe described by the concor- 
dance model of cosmology. Massive galaxies at high z become more and more important to 
study. In fact, the availability of multi- wavelength data from new generations of deep surveys, 



dates up to z ~ 4 — 6.5 


(Giavalisco et al. 200^ 


4: iMobasher et al. 


2005: McLure et al. 


2006 




Yan e 


t all 2006: 


Rodiehiero et allboOT: Bouwens et al.lboOT: Wiklind et al 


I2OO8: Mancini et al. 


2009; 


Stark et a 


. 2009: ] 


Uobas 


ler & Wiklind 


2010l: ] 


Dahlen et al 


2010: ( 


Capak et al. 


2011 


) 


or even z 1{] 


Bouwens et al. 


2011). However, still verv few is known about the phvsical 



processes underlying the (putative) presence of such massive systems at these high redshifts. 
In fact, the mere existence of them could seem pa radoxical with i n a direct interpreta tion of 



the hierarchical structure formation scenario (e.g. iToomrd 119771 : IWhite fc Reeslll978l ). Fur- 
ther, such possibl e contradictions are not necessarily al l eviated by th e competing monolithic 
collapse scenario (jEggen. Lynden-Bell fc Sandagelll962l : lLarsonlll974j ) and hence the question 



about the existence (and the properties) of massive galaxies at high-redshift remains open. 

Recently, a scenario has emerged to explain massive galaxy formation that shares charac- 
terisctics of both t he aforementioned classica l scenarios, but is nevertheless di fferent. Indeed, 



analytical models (ISalvador-Sole et al.ll2005l ). as well as N-body simulations ( IWechsler et al. 
2002; Zhao et al. 20031 ). have shown that two different phases can be distinguished along halo 
mass assembly: i) first, a violent, fast phase, with high mass aggregation (i.e., merger) rates, 
ii) later on, a slow phase, where the mass aggregation rates are much lower. Hydrodynamical 
simulations have confirm ed this scenario and its imp li cation s for prop erties of massive g alac- 
tic objects at low z, see iDommguez-Tenreiro et al.l ( 120061 ) , see also lOser et al.l ( I2OIOI ) and 
Cook. Lapi fc Granatol J2009h . Concerning high ^s, it has been shown that the fast phase 
has the characteristics of a multiclump collapse, where mergers involve very low relative an- 
gular momentum, and, in fact, the y are induced by the collapse of flow convergence regions 



displaying a web- like morphology (IDommguez-Tenreiro et al.l 120101 ). 



In this Letter we investigate the high mass end of galactic stellar objects at high-redshift 
(z = 4, 5 and 6) obtained by means of self-consistent cosmological simulations within a 
volume large enough to account for the proper treatment of the large-scale structure yet 
simultaneously capturing all the relevant small-scale (baryonic) physics. Not only do we 
investigate their mere presence, we also study whether they had enough time to dynamically 
relax at such high redshifts. To this end, we have focused on the intrinsic mass (as opposed 
to luminosity) as well as structural and kinematical properties of these objects at their halo 
(i.e. virial radius) and stellar /galactic scale (to be defined below) as fingerprints of the 
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physical processes involved in their assemblyll] Specifically, we investigate the appearence of 
samples of high-z massive galactic objects with dynamical planes tilted rela tive to the virial 



plane, and link their underly ing formation physics to the Adhesion Model ( iGurbatov et al. 
19891 : IVergassola et alJEooi ). 



Our results here are an extension to high er zs of previous studies on the fundamental 
plane (FP) at 2; = (lOnorbe et al.ll2005U2006l ). In these two papers the different possibilities 
causing the tilt of the FP relative to the virial plane are analyzed in detail. It is shown that 
that if both the viri al mass to luminosity ratio, M^ir/L, and the mass structure coefficient cJJ[ 
(see Eqs. 2 and 4 in lOnorbe et al.l ( 120051 )) are independent of mass, then no FP tilt would be 
measured. Here, because mass is considered instead of luminosity and no projection effects 
are taken into account, we have instead analyzed the mass dependence of the ratio M^^^^/M^ij. 

where M^^^'^ and af!'j^'^ are the 



and the mass structure coefficient 



e.bo/' 



bo 



stellar mass and the 3D stellar velocity dispersion at the galactic scale, respectively, and 



^star 
e.bo 



is the 3D stellar mass effective radius also defined at the galactic scalqj. We need to 
stress that a mass dependence of either of these quantities automatically implies a tilt of the 
dynamical plane relative to the virial plane. 



2. Structural and kinematical Properties of massive Objects at high-2; 



The simulations used here are part of the GALFOBS project. They are iV-body + SPH 
simulations that have been performed usi ng an OpenMP parallel version of the DEVA code 
( ISerna. Dommguez-Tenreiro. fc SaizI 20031) and the methods for star formation and cooling 
described in iMartmez-Serrano et al.l ( l2008l ). The DEVA code pays particular attention to 
ensure that conservation laws (e.g. momentum, energy, angular momentum, and entropy) 
hold as accurately as possibl^. Star formation is implemented through a Kennicutt-Schmidt- 
like law with a density threshold pthres and a star formation efficiency of c. The values of 
these parameters implicitly account for star formation regulation by discrete energy injection 
processes. 



^ We stick to 3D properties as w ell as mass (instead of luminosity) for two reasons: i) projection effects 
add noise in the statistical analysis (|Oiiorbe et al.ll2006l ) and ii) we are not aiming at providing observables 
but rather at understanding the physical processes involved in the formation of these objects. 



^Please refer to Table 1 of lOnorbe et al.l (|2006l ) where our nomenclature and definitions are more thor- 
oughly introduced. 

■^This in particular implies that a double loop in the neighbour searching algorithm must be used, which 
considerably increases the CPU time. 
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The main simulation was carried out in a periodic cube of 80 Mpc side length using 512^ 
baryonic and 512^ dark matter particles with a gravitational softening of eg = 2.3 kpc and 
a minimum hydrodynamical smoothing length half this value. The cosmology applied was a 
ACDM model whose parameters as well as those of the field of primordi al density fluctuation s 
(i.e., initial spectrum) have been taken from CMB anisotropy data0 (jPunkley et al. 2009 ). 
with Qm = 0.295, Qb = 0.0476, Qa = 0.705, h = 0.694, an initial power-law index n = 1, 
and (Tg = 0.852. The mass resolution is rrihar = 2.42 x lO^M© and rridm = 1-26 x IO^Mq and 



the star formation parameters used were pthres = 4.79 x 10 



gem 



and c = 0.3. 



When analyzing galaxy formation in numerical simulations it is desirable to verify that 
the objects in the simulation are consistent with observations at low z's. Due to the extreme 
CPU consumption by hydrodynamical forces, this is not yet possible for the main GALFOBS 
simulation. As a way out, we ran three sub- volumes of the main cube using a "zoom 
approach". In this approach the gravitational forces have been calculated for the full box 
whereas the hydrodynamical forces (which are exclusively local) were only computed in a 
sub-box of side length 26 Mpc. These three sub-volumes have been analysed at redshift 
z = showing that we indeed obtain galaxy populations i n agr eement with low-redshift 
observations, as we had also previously shown in lSaiz et al. I ( 120041 ) using the same approach 
yet smaller simulation boxes. 

^Halos in our simula tions are identified by the OpenMP+MPI halo finder AHllfl ( KnoUmann Sz Knebe 



20091 ) as well as SKID (IWeinberg et al.lll997l ). and their respective results have been cross- 
compared to check for com pleteness. The halo scale of these objects is defined by the virial 
radius (rvir) based upon the iBryan fc NormanI (119981 ) fitting function to determine the over- 
density threshold. The so-called galactic scale has been based upon material (stars) inside 



a sphere of radius 



0.15 X Tvir, a scale separating the baryon from the dark matter 



domination (IBailin et al.ll2005l ). This automated procedure has been tested by comparing 
with individually determined limiting stellar sizes of several hundred of objects based upon 
their 3D visualization as well as their 3D stellar density profiles. We further asked that our 
objects are not involved in violent events, either at the halo or at the galactic scale. To 
exclude this kind of objects at the halo scale, we have used the form factor, cp, defined via 



the virial relation cp = (GMvir)/([o" 



totl2^tot 



e,h 



), where is the velocity dispersion and r 



tot 
e,h 



the half-mass radius at the halo scale, and we asked it to be within the expected interval 
(1.9, 2.5) for virialized objects (IBinney fc Tremaind 120081 ) and of the order of unity if we 



4 



http: / /lambda. gsfc.nasa.gov/product /map/current /params /lcdin^z_lens_run_wmap 



5 _bao_snalLlyapost . cfm 

'^AHF can be freely downloaded from http://popia.ft.uam.es/AMIGA 
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use Tvir instead of r^"^. The same procedure has been employed on the galactic object scale 
using - in analogy - the parameter cf^^ = (GM^Jf )/([cr^Q^'']^rf{^o)- Again, objects outside 
a certain range (based upon a manually gauged subsample of 200 objects for each z) have 
been discarded. Putting a mass threshold of M^^^^ > 10^^ Mq, our final samples consists of 
137, 521 and 1315 galaxies aX z = 6, z = 5 and z = 4, respectively, not involved in violent 
events at any scale. 

Our first result is in fact the mere existence of these samples of high-redshift relaxed 
massive galaxies. To understand their origin, we first quantify the correlation and inter- 
relation, respectively, between their mass (Mvir and M^^^'^), size (half-mass radii rg°j^ and 
rf^l) and velocity dispersior^ (crj;°* and (J^^'^) both at the halo and at the stellar scale 
using the following variables: Ehaio = logigM^ir, rhaio = \og^orl°^, vuaio = \og^Q(J^^\ and 
Estar = logioM^*''^ r.tar = logiQr^^bo' ^star = hg^Q a^^^"" . Wc list the average values E, f, and 
V in Table [1] where we can observe a mild increase of Vhaio alongside a decrease of Vhaio while 
Ehaio remains constant: as the Universe expands, the objects become on av erage less and 



less compact due to the decrease of the global density ( iPadmanabharu 119931 ) . We basically 
observe the same phenomenon on the stellar scale, however, accompanied by a moderate 
increase in Egtar- In fact, the ratios rl°^/rf'^^l and crl^^/a^^^ show an scaling behaviour as a 
function of either M^^j. or M^^^^ (see Table 2 and below). 

Going one step further, we search for planes in the [E, r, v) space by performing a 
princip al component analys is (PGA) of all samples. It is made in 3D to circumvent projection 



effects (lOhorbe et al.ll2006[ ). We have found that at all redshifts one of the eigenvalues of the 
PGA is considerably smaller than the others, so that (massive) objects populate a flattened 
ellipsoid close to two-dimensional, both at the halo scale and at the stellar object scale: 

Es-E, = «f (r, - f,) + 7f (t;, - v,), (1) 

where s refers to the scale of the object, i.e. halo or star. Tabled] also lists the values of the 
parameters and of these planes as well as their bootstrapping errors. We find that 
at the halo scale the planes are close to the virial plane (VP, defined by (cr{;°*)^ = GM^i^/r^^^ 
^haio — — 2), as expected for well-defined haloes. At the stellar scale we also 

find planes to which we refer as the intrinsic dynamical planes (IDPs) and whose observed 
manifestation is the Fundamental Plane. 

To better view the IDP's, their relation to the VP and any possible evolution we plot 
them in Figure [1] for z = Q (green), 2; = 5 (red), and z = A (blue). Points are the actual 
data for all the massive objects, shown in a projection where the z = A data are edge-on. 



^We stress that all our objects here are velocity dispersion supported 
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Ellipses stand for the corresponding projections of the 1-sigma 3D ellipsoids (full lines) or 
3-sigma (blue dashed line). The centers of the ellipses are the corresponding projections of 
the ellipsoid centers. Straight lines have the same directions as the major axes of the ellipses 
resulting from projections of the VPs ellipsoids. Two important results arise from this plot 
and Table [H First, high-z massive galaxies are on IDP's which are clearly tilted relative to 
the VP. Second, we observe a mild evolution of the IDP between z = 6 and z = 4, primarily 
driven by changes in the average values Egtar, fstan Vgtar and not the plane parameters a^^^ 
and 7fi^^. 



3. Discussion & Conclusions 

To shed more light onto the tilt of the IDP with respect to the VP, and following 
the discussion in ^ we have first checked if there is a mass homology breaking, that is, 
if the parameter depends on M^^^^. We have calculated its trend with stellar mass 
(logM^*^^'^ oc P^^logc^') and listed the best fit jS^l in Table H Within the error bars the 
correlation is consistent with zero. This means that stars accomodate the product of their 
spatial and velocity dispersion distributions (i.e., rf'^bo'^\^o^) according to Mvir. Second, we 
have checked whether the mass ratio M^^^'^ /M^ij. correlates with M^*^''. And in fact, we find 
that this ratio decreases for increasing stellar mass at any given z (cf. Table [2]). Therefore, we 
ex pect the IDP to be t ilted against the VP (as discussed in the Introduction and explained 



in lOfiorbe et al.l l2005l . 120061 ). We further compared our results against a simulation with 



different star formation parameters and found no difference concerning the IDP tilt. 

But how can we understand these trends with respect to scenarios of galaxy formation? 
In order to answer this question we need to additionally consider the ratio of hot and cold 
baryon mass inside (i.e. M^^/M^^) function of the mass scale. First we note that 
mass assembly of t he ob jects we analyze is dominated by cold accretion mode, in consistency 



with iKeres et al.l ( 120091 ) results with the entropy-conserving GADGET-2 code. Now, the 
best-fit parameters to the scaling relation logM^Jf'' oc P log M^^/M^ are given in Tabled 
again. There we find that the fraction of hot over cold baryons increases very significantly 
as we go to higher masses. We can further acknowledge from Table [2] that the overall baryon 
fraction = Mj'j'^YMvir does not depend on M^Jf . Both these results taken together imply 
that massive haloes have proportionally less cold gas available to be accreted from the halo 
and transformed into stars than less massive ones. This explains the trend of M^^^"^ /M^[^ 
with the mass scale found above. Further, it is worth noting that all the /3 slopes in Table [2] 
change from one z to another only within their errors. 

Now, why is more hot gas relative to cold gas enclosed within the virial radius as Mvir 
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increases? To answer to this question we have to recall how massive galaxies assemble their 
mass. Very briefly, our simulations show that massive galaxies form from gaseous mass 
elements enclosed by overdense subvolu r nes within the simulati on box. As predicted by the 
Adhesion Model (iGurbatov et al.l Il989l : IVergassola et al.l Il994| ) we have found that gas is 
bi-phasic. Indeed, at a given time a distinction can be made between singular gaseous mass 
elements (as those that have already been involved in caustic, i.e., singularity, formation 
at this time) and regular ones (those that have not yet been trapped into a caustic and 
tend to be of low density). We have also found that, from a global point of view, mass 
elements are dynamically organized as a hierarchy oi flow convergence regions (FCRs), that 
is, attraction basins for mass flows. At high z FCRs undergo fast contractive deformations, 
that violently shrink them, transforming most of the cold, densest gaseous mass elements 
they contain into stars and heating the diffuse component there. Due to its low density, 
this component, once heated te nds to keep hot along ev olution, and forms shock fronts 
that expand, in consistency with iBirnboim fc Dekell ( l2003l ) results. We refer the reader to 



Dommguez-Tenreiro et al.l (120101 ) for a more elaborate discussion. 



In the simulations analyzed in this Letter, we have witnessed events occurring along 
the fast phase of massive galaxy formation, see §T] very fast mass assembly, dissipation 
and star formation rates ensuing FCR contractive deformations. These contractions act 
on dynamical timescales that are short because we have high overdensities where massive 
galaxies are about to form, therefore explaining the presence of massive objects in a young 
Universe. Additionally, such violent FCR contractions tend to swallow the mass close to 
them, severely limiting the amount of mass available to be further assembled after they 
occur. This would explain why a fraction of the objects we have identifled are not dynamically 
disturbed. We have also seen in the simulations the gravitational gas heating due to these 
violent dynamical events, that partially transform the ordered mechanical energy involved 
in contractions into thermal energy and pressure. This is a crucial point for understanding 
the tilt of the IDP's at high z^s. To be quantitative, recall for example that a system must 
get rid of an amount of energy equ al to its binding energy as it collapses from inflnity and 
virializes ( iBinney Sz Tremainell2008l ). This binding energy per unit virial mass increases with 
halo mass as M^/^, so that at assembling a galaxy, the more massive it is the more energy per 
unit mass is available to heat and pressurize the gas at the corresponding FCR contraction. 
Otherwise, as explained above, after these violent events most of the heated low density gas 
elements remain hot. This implies that more hot gas relative to cold gas is enclosed within 
r^ir as Mvir increases, as we have found. Therefore, we can conclude that the origin of the 
IDPs tilt relative to VPs lies in that gravitational gas heating processes are more effective 
as the mass of the halo increases and that there is not mass homology breaking. Finally, let 
us stress that the same physical processes act along the fast phase, as the slopes in Table [2] 
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do not change with zs. 

Summing up, the processes involved in high- 2; massive galaxy formation are: FCR con- 
tractions (approximatively equivalent to collapse) acting on a bi-phasic gas, induced by 
singularity formation in terms of the Adhesion Model approximation; the ensuing transfor- 
mation of the ordered mechanical energy of contraction into velocity dispersion, and then 
partially into thermal energy and gas pressure, on the same timescales; and dense gas ele- 
ments shrinkage, cooling and their transformation into stars. Energy injection is unlikely to 
substantially change the processes responsible for this high z FP tilt, because, as explained, 
they have to do with caustic (i.e. singularity) formation. We conclude that the violent 
processes described above are responsible for having: 1) massive objects at high redhisft, 2) 
hot gas coronae, 3) less cold gas to form stars as the mass scale increases, because we have 
more gas heated, implying IDP's tilted relative to VP's, among other results. We see that 
the same processes are responsible for obtaining massive stellar objects shortly after the Big 
Bang as well as having them lying on IDPs: fast FCR contractions at different scales are the 
engine driving them. 

We thank the anonymous referee for useful suggestions that improved the manuscript. 
We thankfully acknowledge to D. Vicente and J. Naranjo for the assistence and technical 
expertise provided at the Barcelona Supercomputing Centre, as well as the computer re- 
sources provided by BSC/RES (Spain). We thank DEISA Extreme Computing Initiative 
(DECI) for the CPU time allowed to CALFOBS project. The Centro de Computacion 
Cientifica (UAM, Spain) has also provided computing facilities. This work was partially 
supported by the DGES (Spain) through the grants AYA2009-12792-C03-02/-03, as well as 
by the ASTROMADRID network (CAM S2009/ESP-1496). AK and JO were supported by 
the MICINN, Spain, through the Ramon y Cajal programme and the " Supercomputacion y 
e-Ciencia" Consohder-Ingenio 2010 project (CSD2007-0050), respectively. 

REFERENCES 

Bailin, J., et al. 2005, ApJ, 627, L17 

Binney J. & Tremaine S. 2008, Galactic Dynamics, Princeton University Press (Princeton, 
New Jersey) 

Birnboim, Y., & Dekel, A. 2003, MNRAS, 345, 349 

Bouwens, R. J., lUingworth, G. D., Franx, M., & Ford, H. 2007, ApJ, 670, 928 



- 9 - 



Bouwens, R. J., et al. 2011, Nature, 469, 504 

Bryan, G.L. & Norman, M.L. 1998, ApJ, 495, 80 

Capak, P. L., et al. 2011. larXiv:1101.3586l 

Cook M., Lapi A., Granato G. L., 2009, MNRAS, 397, 534 

Dahlen, T., et al. 2010, ApJ, 724, 425 

Dommguez-Tenreiro R., Ofiorbe J., Saiz A., Artal H., & Serna A. 2006, ApJ, 636, L77 

Dommguez-Tenreiro R., Ofiorbe J., Martmez-Serrano F. J. & Serna A. 2010 MNRAS, ac- 
cepted 

Dunkley, J., et al. 2009, ApJS, 180, 306 

Eggen O. J., Lynden-Bell D., Sandage A. R., 1962, ApJ, 136, 748 
& Dave R. 2001, ApJ, 562, 605 

Giavalisco M., Dickinson M., Ferguson H. C., et al. 2004, ApJ, 600, L103 

Gurbatov S. N., Saichev A. I., Shandarin S. F., 1989, MNRAS, 236, 385 

Keres, D., Katz, N., Fardal, M., Dave, R., & Weinberg, D. H. 2009, MNRAS, 395, 160 

KnoUmann, S. R., & Knebe, A. 2009, ApJS, 182, 608 

Larson R. B., 1974, MNRAS, 166, 585 

Mancini, C., Matute, I., Cimatti, A., Daddi, E., Dickinson, M., Rodighiero, G., Bolzonella, 
M., & Pozzetti, L. 2009, A&A, 500, 705 

Martmez-Serrano, F. J., Serna, A., Domfnguez-Tenreiro, R., & MoUa, M. 2008, MNRAS, 
388, 3 

McLure, R. J., Cirasuolo, M., Dunlop, J. S., et al. 2006, MNRAS, 372, 357 
Mobasher B. et al, 2005, ApJ, 635, 832 

Mobasher B., Wiklind T., 2010, in The Impact of HST on European Astronomy, Astrophysics 
and Space Science Proceedings, F. Duccio Macchetto ed.. Springer, Netherlands 

Ofiorbe, J., Domfnguez-Tenreir,o R., Saiz, A., Serna, A., Artal, H. 2005, ApJ, 632, L57 

Onorbe, J., Domfnguez-Tenreiro, R., Saiz, A., Artal, H., Serna, A. 2006, MNRAS, 373, 503 



-lo- 
oser, L., Ostriker, J. P., Naab, T., Johansson, P. H., & Burkert, A. 2010. larXiv:1010.1381l 

Padmanabhan, T. 1993, Structure Formation in the Universe, by T. Padmanabhan, 
pp. 499. ISBN 0521424860. Cambridge, UK: Cambridge University Press, June 1993. 

Rodighiero, G., Cimatti, A., Franceschini, A., Brusa, M., Fritz, J., & Bolzonella, M. 2007, 
A&A, 470, 21 

Saiz, A., Dommguez-Tenreiro, R., & Serna, A. 2004, ApJ, 601L, 131 
Salvador-Sole E., Manrique A., Solanes J. M. 2005, MNRAS, 358, 901 
Serna, A., Dommguez-Tenreiro, R., & Saiz, A. 2003, ApJ, 597, 878 

Stark, D.P., Ellis, R.S., Bunker, A., Bundy, K., Targett, T., Benson, A., Lacy., 2009, ApJ, 
697, 149 

Toomre A., 1977, in The Evolution of Galaxies and Stellar Populations, eds. B. Tinsley & 
R. Larson (New Have, CN: Yale Univ. Press) 

Vergassola M., Dubrulle B., Frisch U., Noullez A., 1994, A&A, 289, 325 

Wechsler R.H., Bullock J.S., Primack J.R., Kravtsov A.V., Dekel A., 2002, ApJ, 568, 52 

Weinberg, D. H., Hernquist, L., & Katz, N. 1997, ApJ, 477, 8 

White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 

Wiklind, T., Dickinson, M., Ferguson, H. C, et al. 2008, ApJ, 676, 781 

Yan, H., Dickinson, M., Giavahsco, M., et al. 2006, ApJ, 651, 24 

Zhao D.H., Mo H.J., Jing Y.P., Borner G., 2003, MNRAS, 339, 12 



This preprint was prepared with the AAS E^TJrjX macros v5.2. 



- 11 - 




Fig. 1. — The z = 4 (blue), z = 5 (red) and z = Q (green) IDPs seen in a projection where 
the z = 4 data are edge-on. Full line ellipses represent the corresponding 1 x a ellipsoids 
seen in the same projection. The 3 x a ellipse is also plotted for the z = 4 sample (long 
dashed line). The centers of the ellipses are the corresponding projections of the ellipsoid 
centers. Data points for all the massive objects in the samples are also plotted as circles 
using darker versions of their respectively colours. Short dashed lines are the projections of 
the major axes of the VPs ellipsoids. M^^^^ in Mq, r^^^J^ in kpc and a^^^ in km x s~^. 



Table 1. Results of PCA analysis. 



Sample 


No. 


E 


f 


V 


a3D 




CBrv 


Halo 


2 = 6 
2 = 5 
2 = 4 


137 
521 
1315 


11.321±0.016 
11.326±0.010 
11.336±0.007 


1.05,3±0.007 
1.118±0.005 
1.188±0.003 


2.319±0.006 
2.284±0.003 
2.248±0.003 


0.930±0.040 
0.822±0.017 
0.798±0.013 


1.906±0.050 
2.008±0.024 
2.053±0.018 


0.0100±0.0005 
0.0092±0.0003 
0.0105±0.0002 


Stellar 


2 = 6 
2 = 5 
2 = 4 


137 
521 
1315 


10.189±0.012 
10.209±0.007 
10.237±0.005 


0.067±0.007 
0.119±0.004 
0.189±0.003 


2.220±0.006 
2.187±0.004 
2.146±0.003 


-0.142±0.091 
-0.002±0.059 
0.077±0.035 


2.096±0.104 
2.040±0.052 
1.994±0.027 


0.0297±0.0021 
0.0290±0.0010 
0.0257±0.0006 



Note. — Column 2: Number of massive galaxies in the sample. Columns 3, 4 and 5: sample mean values of the E, 
r and v variables (log M©, log kpc and log km x respectively). Columns 6 and 7: coefficients of the IDP plane. 
Column 8: IDP orthogonal scatter in the E,r and v variables. Errors have beend obtained from a bootstrapping 
analysis of the samples. 



Table 2. Direct Fits. 



CF 


-0.071 ±0.009 


-0.076 ±0.014 


-0.045 ± 0.033 


1.0173 


1.0500 


0.7267 


„vir 


0.059 ± 0.032 


0.019 ± 0.066 


0.074 ±0.142 


0.8883 


1.2765 


0.6917 




-0.121 ±0.023 


-0.131 ±0.046 


-0.137 ±0.096 


1.4470 


1.5277 


1.5876 


„tot /„star 
^e,h/'^e,bo 


0.252 ± 0.025 


0.225 ± 0.046 


0.256 ±0.110 


-1.5761 


-1.3013 


-1.6227 




-0.306 ± 0.029 


-0.258 ± 0.057 


-0.222 ±0.116 


-2.0301 


-1.5159 


-1.1264 




0.789 ± 0.022 


0.829 ± 0.040 


0.877 ±0.091 


-8.8637 


-9.2277 


-9.6826 




-0.061 ±0.006 


-0.048 ± 0.009 


-0.032 ± 0.142 


-0.1594 


-0.2892 


-0.4581 



Note. — Correlation between various properties X and M^^^ as derived from fitting log Af^^*"" = 
/9 logX + (p. Errors stand for a 97.5% confidence level intervals. 



